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Abstract 

The navigation algorithms for the Sea- viewing Wide Field-of-view Sensor (SeaWiFS were designed to meet the 
requirement of 1 -pixel accuracy — a standard deviation (<r) of 2. The objective has been to extract the best 
possible accuracy from the spacecraft telemetry and avoid the need for costly manual renavigation or geometric 
rectification. The requirement is addressed by postprocessing of both the Global Positioning System (GPS) 
receiver and Attitude Control System (ACS) data in the spacecraft telemetry stream. The navigation algorithms 
described are separated into four areas: orbit processing, attitude sensor processing, attitude determination, 
and final navigation processing. There has been substantial modification during the mission of the attitude 
determination and attitude sensor processing algorithms. For the former, the basic approach was completely 
changed during the first year of the mission, from a single-frame deterministic method to a Kalman smoother. 
This was done for several reasons: a) to improve the overall accuracy of the attitude determination, particularly 
near the sub-solar point; b) to reduce discontinuities; c) to support the single-ACS-string spacecraft operation 
that was started after the first mission year, which causes gaps in attitude sensor coverage; and d) to handle data 
quality problems (which became evident after launch) in the direct-broadcast data. The changes to the attitude 
sensor processing algorithms primarily involved the development of a model for the Earth horizon height, also 
needed for single-string operation; the incorporation of improved sensor calibration data; and improved data 
quality checking and smoothing to handle the data quality issues. The attitude sensor alignments have also 
been revised multiple times, generally in conjunction with the other changes. The orbit and final navigation 
processing algorithms have remained largely unchanged during the mission, aside from refinements to data 
quality checking. Although further improvements are certainly possible, future evolution of the algorithms is 
expected to be limited to refinements of the methods presented here, and no substantial changes are anticipated. 


1. INTRODUCTION 

The navigation processing for the Sea-viewing Wide 
Field-of-view Sensor (SeaWiFS) data is performed as part 
of the level -0 to -la conversion. The level -0 to 'la soft- 
ware extracts and converts the required telemetry from the 
data stream and passes it to the navigation code, which 
produces per-scan-line spacecraft position and instrument 
pointing information. The output of navigation is stored 
in the level -la data products for use by downstream pro- 
cessing. 

The navigation code is composed of two, largely in- 
dependent subsystems: orbit processing, which filters the 
data from the onboard global positioning system (GPS) 
receiver to produce orbit vectors; and attitude processing, 
which filters the spacecraft attitude control system (ACS) 
telemetry and instrument tilt telemetry to determine the 
SeaWiFS sensor orientation. 

The remainder of this section defines constants, refer- 
ence frames, and transformations which are used in the 
algorithm descriptions. 

1.1 Constants 

The following constants are defined here for later use: 

■ R&, Earth equatorial radius (6,378.137km); 

■ Rm, Earth mean radius (6,371km); 

■ /, dimensionless Earth flattening factor (1/298.257); 

■ uJE) Earth rotation rate (7.2921 1585494 x 10“ 5 s -1 ); 


■ wo, nominal Orbit angular rate (27r/5940); 

■ G m , Earth gravitational constant (398600.5 km 3 
s~ 2 ); and 

■ J 2 , dimensionless Earth gravity field perturbation 
term (1.08263x 10“ 3 ). 

1.2 Reference Frames 

In order to describe the navigation algorithms, several 
basic reference frames (all frames have orthonormal axes) 
are defined below. 

a. Earth- Centered Inertial (ECI): This reference frame 
has its origin at the Earth’s center and is inertially 
fixed. The axes are defined as: r on the equator at 
the vernal equinox; z at the North Pole; y orthogo- 
nal to z and x in the right-hand sense. 

b. Earth-Centered Earth-Fixed (ECEF): This refer- 
ence frame also has its origin at the Earth’s cen- 
ter and rotates with the Earth. The axes are de- 
fined as: x at 0° latitude and longitude (Greenwich 
meridian at the equator); y at 0° latitude and 90° 
longitude; and z at the North Pole (also known as 
Earth-Centered Rotating, or ECR). 

c. Orbital: This frame has its origin at the spacecraft 
position and is defined as: r-axis along the geodetic 
nadir vector; y-axis perpendicular to x and opposite 
the spacecraft velocity vector; and z-axis toward the 
orbit normal. 


1 



Navigation Algorithms for the SeaWiFS Mission 


d. Spacecraft: This frame has its origin at the space- 
craft position and is defined in the same sense as 
the orbital frame; during normal spacecraft opera- 
tions the two are related by a small rotation which 
is specified by a pitch (#), roll (<£), and yaw (ip) 
Euler angle sequence. 

e. SeaWiFS Base: This frame is nearly the same as 
the spacecraft frame except for a small, constant 
misalignment transformation. 

f. SeaWiFS Instrument: This frame is related to the 
SeaWiFS base frame by the tilt angle, i.e., a ro- 
tation about the 2 -axis, which is common to both 
frames. 


O y = 0 z X 0 Xi (5) 

where Vc is the orbit velocity vector corrected for the 
Earth’s rotation rate, 


Vc = 


V x -U> E Vy 
V y +w E V x . 

v z 


The full transformation is given by 


O = 


\on 

0 T y ■ 

loll 


( 6 ) 


(7) 


1.3 Transformations 


where O is the ECEF-to-orbital transformation matrix. 


The transformations among the ECI, ECEF, orbital, 
and spacecraft reference frames defined above are described 
in the following sections. 

1.3.1 ECI-to-ECEF Transformation 

The ECI-to-ECEF transformation consists of a rotation 
of the Greenwich hour angle (£) about the z-axis (Patt and 
Gregg 1994): 


E 


cos£ 
— sin£ 


0 


sin£ 0 
cos £ 0 . 

0 1 


( 1 ) 


1.3.2 ECEF-to-Orbital Transformation 

The ECEF-to-orbital transformation js derived from 
the orbit position and velocity vectors, P and V t in the 
ECEF frame (Patt and Gregg 1994). First, the nadir vec- 
tor at the spacecraft position is calculated. This is not a 
closed-form calculation, but a very good approximation is 
given by (Patt and Gregg 1994): 


O r = 


(l-/p) 2 ^xl 

(Wp) 2 ^, 

Pz 


1 

v /p2 + (l-/p)^(P I 2 + P y 2 )’ 


( 2 ) 


where 6 X is the geodetic nadir vector, and fp is the effec- 
tive flattening factor at the spacecraft position. The latter 
is computed as follows: 

(i~fp ) 2 = + \Eiz. ?±L } ( 3 ) 


where Rm is the Earth mean radius. The remaining axes 
of the orbital frame are defined as: 


and 


Oz 


Vc x O x 

|Vc x Ox 


(4) 


1.3.3 Orbital-to-Spacecraft Transformation 


The orbital-to-spacecraft transformation is defined by 
an Euler angle sequence. The order of rotations is pitch, 
roll, and yaw, which corresponds to z, y, and x in the 
spacecraft or orbital frame. The transformation is given 

by 

B = (8) 


where 


and 


A^ 

A# 


1 0 0 

0 cos ip sin ip , 

0 — sin ip cos ip 

cos <P 0 sin <p 

0 10, 

— sin (p 0 cos <p 


Aq 


cos 0 
sin# 


0 


— sin 9 0 
cos 9 0 . 

0 1 


(9) 

( 10 ) 


( 11 ) 


2. ORBIT PROCESSING 

The orbit processing algorithms use the telemetry from 
the GPS receiver onboard OrbView-2f (OV2) to compute 
the spacecraft position at the time of each SeaWiFS scan 
line. The approach described here was designed to satisfy 
the following objectives: eliminate excursions and outliers 
in the data; and span intervals of missing, degraded, or 
corrupted data. 

The approach is to use the GPS vectors from the teleme- 
try as observations to be fit to a high-fidelity orbit model. 
The OV2 GPS receiver telemetry includes orbit position 
and velocity vectors in the ECEF frame, associated time 
tags, and the number of GPS satellites tracked. The telem- 
etry is updated every 10 s. 


f The spacecraft carrying the SeaWiFS instrument was origi- 
nally named “SeaStar” ; it was designated “OrbView-2” after 
launch. 
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The orbit model is based on the Artificial Satellite 
Analysis Program (ASAP), which was developed at the 
Jet Propulsion Laboratory (Kwok 1987) and is now avail- 
able through the NASA Computer Software Management 
and Information Center (COSMIC) library of public do- 
main software. ASAP is a Cowel orbit integration pro- 
gram which includes an Earth gravity model of up to 36 in 
degree and order, solar and lunar attractions, and atmos- 
pheric drag. 

The processing for each interval is initialized using the 
results from the previous interval, which are saved as mean 
orbital elements at the ascending equator crossing. The 
orbit is integrated through the data interval using ASAP, 
and the vectors are extracted at the times of the GPS 
data samples. The position differences and partial deriva- 
tives, with respect to the initial elements, are computed 
and used to update the elements via a least-squares esti- 
mation. These calculations are iterated if necessary until 
the solution converges. The final elements are then used 
to integrate the orbit through the data interval, and the 
orbit vectors are interpolated to the SeaWiFS scan line 
times. The details of the orbit processing are presented in 
the following sections. 

2.1 Orbit Integration and Interpolation 

The ASAP orbit model and integration algorithm are 
described in Kwok (1987) and will not be repeated here. 
The model is initialized using mean elements, which con- 
sist of the semimajor axis (a), eccentricity (e), inclination 
(z), right ascension of the ascending node (/), argument of 
perigee (u>), and mean anomaly (m). A 60 s output interval 
is selected. Additional force terms (e.g., atmospheric drag, 
other-body attractions) are based on input parameters and 
models. 

The vectors are generated by ASAP in the ECI frame 
and are transformed to the ECEF frame using the ECI-to- 
ECEF transformation from Sect. 1.2. The vectors are in- 
terpolated to intermediate times using cubic interpolation, 
which combines the position and velocity of each compo- 
nent as the displacement and first derivative of a scalar 
function. If the time tags of two successive vectors are 
arbitrarily assigned values of 0 and 1 , then the cubic poly- 
nomial coefficients are computed as follows: 


Cb = Pi, (12) 

ci = ^ At, (13) 

c 2 = 3p 2 - 3 Pi - 2 V\ At - V 2 A t, (14) 

and 

C 3 = 2 Pi — 2 P 2 F Vi At F V^At, (15) 


where Pi and P 2 are the position vectors at times t\ and 
* 2 ; y\ and V 2 are the velocity vectors; At = t\ — 12 is the 
time difference (nominally 60 s) ; and cq , ci , C 2 , and C 3 are 


the vectors of cubic polynomial coefficients. The vectors 
at intermediate times are then calculated as: 


P = Cq + ci£ F c 2 t 2 F c^t 3 , 
Ci F 2 c2^ F 3 C 3 P 


and 


f = 


At 


(16) 

(17) 


where P is the position vector at time £ 5 ; similarly for 
V\ and £ = (£5 - ti)/A£ is the relative. time, for which 
£5 is the desired sample time. This algorithm is used for 
interpolating the vectors to both the GPS sample times 
and the SeaWiFS scan line times. 

2.2 Partial Derivatives 

The least-squares algorithm requires that partial deriva- 
tives of the orbit position components be computed at each 
data point with respect to changes in the initial mean ele- 
ments. Calculation of analytic partial derivatives account- 
ing for all force model terms would be complex. For small 
changes to the initial elements, however, most terms can 
be neglected, and good results are achieved by considering 
only the central body force and the J 2 (oblateness) term. 
The partial derivatives can be expressed in terms of the 
initial elements, the position and velocity vectors, and the 
elapsed time from the element epoch. The derivatives with 
respect to each element are computed as. follows (note that 
in these equations all angles have units of radians). 

2.2.1 Semimajor Axis 

The partial derivatives with respect to the semimajor 
axis are computed as: 

dP _ P - \{t - tp)V 


da 


\p\ 


(18) 


where P is the orbit position vector, V is the orbit velocity 
vector, £ is the time tag for the vector and £0 is the initial 
element epoch. 

2.2.2 Eccentricity 

The partial derivatives with respect to the eccentricity 
are computed as: 


dP 

He 


2\P\ V sin v 

l?l 


— Pcosu, 


where z/, the true anomaly, is computed as 

dw 


v = o — w 


{t - t Q ) 


dt 


(19) 


( 20 ) 


The time of the orbit vector is £; to is the epoch of the 
orbital elements; and o is the orbit angle measured from 
the ascending node: 


o = tan’ 


1 sin(t 


)(P x cos(J) + P v sin(i)) 


( 21 ) 
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using the signs of the numerator and denominator to eval- 
uate the arctangent over the range of — n to 7r; and 

^ | J 2 y/G^R 2 E a- 7/2 (2 - ^sin 2 (i)), (22) 


where 


L = 


PxV 

\Px9\’ 


is the orbit plane normal vector; 


(28) 


is the time derivative of the argument of perigee. 


? = | h VG^R 2 e a_7/2 sin(t) (t - t 0 ), (29) 

oi 2 


2.2.3 Right Ascension of the Ascending Node 

The partial derivatives with respect to the right ascen- 
sion of the ascending node are computed as: 

^ t x P, (23) 

dl 

where 2 is the geocentric North Pole vector (0,0, 1) T . 


2.2.4 Argument of Perigee 

The partial derivatives with respect to the argument of 
perigee are computed as: 


or ^ . 

— — = P sm v — 
ow 


2\P\V cosv 

e W\ 


(24) 


Note that for near-circular orbits, the total effect of a 
change in the argument of perigee is nearly the same as 
for the mean anomaly, with the difference being a small 
periodic term. To avoid degeneracy, this equation includes 
only the periodic term, and the final value of the mean 
anomaly is corrected to account for the computed change 
in the argument of perigee. 


is the derivative of the right ascension of the ascending 
node with respect to i\ and 

= 12J 2 VG^Re « -7/2 sin(i) ( t - t 0 ), (30) 

Ol 

is the derivative of the mean anomaly with respect to i . 

These equations are complete as given; however, some 
computational efficiency is gained by taking advantage of 
relationships between orbital elements and the orbit vec- 
tors themselves, thereby reducing the number of trigono- 
metric function evaluations. Specifically, 


sin(Z) = 

Lx 



(31) 

+ 

Vy 


cos(Z) — 

-L, 

/ 



(32) 

+ 

% 


sin(z) = 

yJH + 

% 

(33) 

cos(i) = 

L„ 


(34) 


2.3 Least-Squares Update Procedure 


2.2.5 Mean Anomaly 

The partial derivatives with respect to the mean anom- 
aly are computed as: 



is the mean motion. 


2.2.6 Inclination 

The partial derivatives with respect to the inclination 
are computed as: 


dP 

di 


r , dl OP 
Lsm{0) + dl ~dl + 


dm dP 
di dm ’ 


(27) 


The initial orbital elements are updated to fit the GPS 
orbit position vectors using a batch- weighted, least-squares 
algorithm, based on the Gauss-Newton procedure (Fallon 
and Rigterink 1978). The 6-element state vector, X, is 
the initial set of orbital elements; the measurement vec- 
tor consists of the individual components of the GPS orbit 
position vectors, with each component treated as an inde- 
pendent measurement; and the observation model vector 
is the corresponding set of position components generated 
from the orbital elements. 

The partial derivative matrix is generated using (18)- 
(34) for each orbit position vector, with the following addi- 
tional calculations performed. First, the derivatives with 
respect to angles (i, Z, u?, and m) are scaled from radians 
to degrees, to be consistent with the units of the elements 
themselves. Second, the derivatives with respect to ec- 
centricity are about two orders of magnitude larger, and 
those with respect to the argument of perigee are about 
the same factor smaller, than those with respect to the 
other elements;, to avoid numerical errors in the calcula- 
tions, these derivatives are scaled by factors of 0.01 and 
100, respectively, with these scalings later removed from 
the updates. 
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Table 1. The conversion factors for attitude sensor measurements. 


Sensor 

Scale Factor 

Offset 

Range 

Sun Sensor 
Horizon Scanner 

2.0022498 x 10~ 4 
5.4931641 x lO" 3 

-2.050304 

0.0 

±2.050304 

0-360° 


Prior to the actual fitting procedure, the input GPS 
vectors are validated. First, the total number of vectors is 
required to be above a specified minimum value; if fewer 
vectors are available, the initial elements are integrated 
through the data interval and used directly for naviga- 
tion, with no update performed. Second, the number of 
GPS satellites tracked for each measured vector must be 
at least four, the minimum value required to obtain a valid 
GPS measurement. Third, the magnitude of the orbit po- 
sition vector is required to be within a specified range. Fi- 
nally, rejection of outliers is performed by computing the 
interquartile range (IQR) for the initial set of measure- 
ment residuals; any vectors with residuals greater than a 
specified multiple of the IQR are rejected. 

The corrections to the state vector are computed as: 

AJ? = [W ± G T G]- 1 [G r (f - ( 5 )], (35) 

where AJ? is the 6-element state update vector; W is 
the 6x6 state weight matrix; G is the matrix of partial 
derivatives of G with respect to ; V is the vector of 
GPS orbit position measurements; and G is the vector of 
corresponding orbit position components generated from 
the initial orbital elements. The state weights are chosen 
to avoid anomalous updates due to reduced visibility on 
some elements for short data intervals [i.e., high-resolution 
picture transmission (HRPT) scenes of 10-15 m], while 
allowing rapid convergence for longer intervals (onboard 
recorder downlinks of approximately 12 h). The measure- 
ment weights are implicitly set to unity. 

Prior to their application, the updates to the eccentric- 
ity and the argument of perigee are scaled by the same 
factors as the associated derivatives (0.01 and 100, respec- 
tively); and the update for the argument of perigee is sub- 
tracted from that for the mean anomaly, as described above 
in the calculation of the derivative with respect to the ar- 
gument of perigee. 

The state vector updates are tested for convergence by 
comparing them with specified tolerances. The tolerances 
are chosen to represent a maximum position change of ap- 
proximately 100m over 12 h: 0.002 km for a; 0.00001 for e; 
0.1 for w ; and 0.001 for i, Z, and m. If the update to any 
element exceeds the tolerance, the entire process (orbit in- 
tegration and interpolation, partial derivative calculation 
and least-squares update) is iterated using the updated 
orbital elements. 


2.4 Final Orbit Processing 

The converged orbital elements are used to generate the 
orbit vectors for use in attitude processing and geolocation 
of the scan lines in the interval. The orbit is integrated 
through the data interval using the ASAP model and the 
converged elements, and orbit vectors are stored every 60 s. 
During attitude processing for individual scenes — global 
area coverage (GAC) orbits, local area coverage (LAC) 
segments, or HRPT scenes — the orbit is interpolated to 
the scan line times using (12)-(17). 

The results of GPS orbit fitting are stored for use in 
future reprocessing, to avoid the need to repeat the fitting 
process, and to initialize the fitting process for the next 
data interval. This is performed by saving the mean el- 
ements at each ascending node crossing in a permanent 
disk file. The converged orbital elements at the start of 
the interval are saved, and the orbit vector at each node 
crossing (i.e., the first vector north of the equator in each 
orbit) is converted to mean elements, which are appended 
to the file. The algorithm for converting the orbit vectors 
to mean elements is included in the ASAP program. 

3. ACS SENSOR PROCESSING 

The OV2 attitude is computed using the spacecraft at- 
titude sensor data, along with the orbit vectors produced 
by filtering the GPS data. The first step is filtering and 
converting the attitude sensor data. The OV2 attitude 
sensor complement consists of three two-axis digital sun 
sensors (DSS) and two horizon scanners (HS). The sample 
rate is once per 10 s in GAC data and once per 2 s in LAC 
and HRPT data. Note that because the ACS telemetry 
are included twice per second in LAC and HRPT data, 
the samples are repeated four times each. The conver- 
sions from digital telemetry to analog measurements for 
the sensors are given in Table 1. The processing converts 
each sensor output to unit vectors in the spacecraft refer- 
ence frame. The output and processing for each sensor are 
described separately. 

3.1 Sun Sensor Processing 

The following sections describe the DSS, their output, 
and the sensor data processing. 

3.1.1 Sensor Description 

The OV2 DSSs each measure the sun’s direction in two 
axes, with a field-of-view (FOV) of ±64° along each axis. 
The three sensors are mounted as follows: DSS- A is on the 
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Fig. 1 . The OrbView-2 digital sun sensor coverage. 


forward side (— y axis) of the spacecraft, with its boresight 
angled 26° upward (toward the — x axis); DSS-B is on the 
rearward (+y) side, also with its boresight angled 26° up- 
ward; DSS-C is on the spacecraft top deck, with its bore- 
sight along the zenith (— x) direction. The DSSs, there- 
fore, provide continuous and overlapping coverage during 
the orbit day. The coverage of the DSSs is illustrated in 
Fig. 1. . 

The output of each sun sensor are the tangents of the 
two projected angles (Chen and Lerner 1978), time tag, 
and status words. The tangents are digitized with a range 
of 0-20,479 counts and are converted to engineering units 
as given in Table 1. The time tags are given in milliseconds 
of the week, and are converted to a time offset from the 
minor frame time tag. The status words include the sun 
presence bit, which indicates whether the sun is within the 
sensor FOV. 

3.1.2 Sensor Data Processing 

The steps in the processing of the sun sensor data are: 
validation, smoothing and interpolation, transformation, 
and averaging. The first two steps are performed indepen- 
dently for each output of each DSS, the two output values 
for each DSS are combined during the transformation, and 
the final stage produces a single vector from all sensors. 

The validation is performed in three steps. First, the 
valid arc for each sensor is determined as the interval dur- 
ing which the sun presence is indicated by the status word. 


Second, the data along each axis are limit-checked against 
an absolute range. Third, the consecutive samples for each 
axis are checked for consistency via a delta limit check. 
Samples which fail any test are rejected from further pro- 
cessing. The limits for the last two checks are adjustable 
parameters. 

The smoothing and interpolation is performed as a sin- 
gle step, to reduce any remaining noise spikes and to inter- 
polate the sensor samples to the SeaWiFS scan line times. 
This is performed using a moving-arc cubic polynomial. 
The polynomial is fit to a fixed number of samples, start- 
ing with the first valid sample, and evaluated at the scan 
line times. The sample range is then shifted by half of 
the number of samples, and the fitting-evaluation process 
is repeated. In the overlapping range, the final output is 
computed as a weighted average of the two polynomials, 
with the weight increasing linearly from the end to the 
center of the sample range for each fit. This process is 
repeated until the range of valid samples has been fitted. 

The transformation converts the output of each DSS 
into unit vectors in the spacecraft frame. The output of 
the sensor is corrected using a linear scale factor and bias: 

tana' = A a tana H- £ a , (36) 

and 

tan/?' = A^tan/? + Bp, (37) 

where tan a and tan /? are the sensor output values, which 
have been validated and smoothed; tana' and tan/?' are 
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the calibrated values; A a and Ap are the scale factors; and 
B a and Bp are the biases for the two axes. The results 
are converted to unit vectors in the DSS frame as follows 
(Chen and Lerner 1978; note that here the sensor x-axis is 
the boresight): 


S D = 


1 

tan a! 
tan/3' 


1 

y/l T tan 2 a ' + tan 2 /J' 


(38) 


nominal values for the phase measurements are about 225° 
and 135° for HS-A and HS-B, respectively, and about 120° 
for the chord measurements for both scanners. 

The output of each horizon scanner consists of the 
chord and phase angles, time tag, and status words. The 
angles are digitized with a range of 0-65,535 counts and 
are converted to engineering units as given in Table 1 . The 
time tags are given in milliseconds of the week, and are 
converted to a time offset from the minor frame time tag. 


and then transformed to the spacecraft frame: 

S s = C S D , (39) 

where C is the DSS-to-spacecraft transformation matrix. 

The final step is to compute the average observed sun 
vector. For each scan line which had sun presence for two 
DSSs, the weighted average of the vectors from the DSSs is 
computed. The weight for each DSS is unity at the center 
of the FOV and decreases linearly to zero at the edges, 
based on the assumption that each DSS is most accurate 
near the center of the FOV. 

3.2 Horizon Scanner Processing 

The following sections describe the HSs, their output, 
and the sensor data processing. 

3.2.1 Sensor Description 

The OV2 HSs measure the Earth’s horizon height (or, 
more precisely, the height of the CO2 layer in the atmos- 
phere) in the scanner reference frame. Each scanner has 
an optical element which scans a conical path (with a half- 
cone angle of 45°) as the sensor rotates. The scanner de- 
termines the rotation angle at which the optical element 
crosses the horizon from space to Earth (in-crossing) and 
vice versa (out-crossing). The scanner determines the ro- 
tation angle between the two crossings (Earth width, or 
chord) and the mid-point relative to a zero- reference angle 
(phase). These measurements are used to determine the 
unit nadir vector in the spacecraft frame. 

The two scanners are mounted on opposite sides of the 
spacecraft lower deck, with their rotation axes toward the 
-z (HS-A) and +z (HS-B) axes. The rotation axes are 
canted 5° toward the +x (downward) axis. Each scanner 
rotates in the right-hand sense about its outward-pointing 
axis. The HS geometry is illustrated in Fig. 2. 

Each scanner has a blanking region , which is a 90° 
range of rotation for which the sensor view is obscured by 
the mounting hardware. The zero-reference for the phase 
measurement is centered within the blanking region, so the 
effective range for horizon crossing measurements is from 
45-315°. The OV2 scanners are mounted with the zero- 
reference rotated 45° from the spacecraft upward (x-axis), 
to enable the scanners to be used both at nadir pointing 
and with the spacecraft pitched 90° for orbit raising. The 


3.2.2 Sensor Data Processing 

The processing of the horizon scanner data is more 
complicated than that for the sun sensors. The reason is 
that the quantity being measured — the height of the CO2 
layer in the atmosphere — is dynamic and varies with orbit 
position, as well as having seasonal and unmodeled varia- 
tions. The total variation from the average is of the order 
of tenths of a degree, which is significant compared to the 
SeaWiFS pixel angle (0.091°), and therefore, the trigger- 
ing height model is critical to the overall accuracy of the 
calculation. 

The initial processing of the scanner data is very similar 
to that for the sun sensors. In the CO2 band (14-16 ^m), 
the scanners can view the Earth for the entire orbit, so 
there is no need to determine the valid data arc. The 
individual measurements are checked using both absolute 
and delta limits, smoothed and interpolated to the scan 
line times, using the moving-arc cubic polynomial method. 

The conversion of the angles to the geocentric nadir in 
the spacecraft reference frame involves several steps: con- 
verting angles to horizon crossing vectors; calculating the 
model horizon angles at the crossings; and solving for the 
nadir vector. 

3.2,2, 1 Horizon Vectors 

The phase and chord measurements are converted to 
incrossing and outcrossing angles, as follows: 

(40) 

and n 

$0 = $s + ^ » (41) 

where $5 is the phase angle, $7 is the chord angle, and 4>/ 
and $0 are the incrossing and outcrossing angles, respec- 
tively, relative to the zero reference for the scanner. 

The crossing angles are calibrated using a look-up ta- 
ble generated from calibration data provided by the sensor 
vendor. A unique calibration table is provided for each 
angle, with measured and actual angles tabulated at 0.5° 
intervals. The actual crossing angles are computed by 
interpolating from the table entries to the input angles. 
A sample calibration table is illustrated in Fig. 3, which 
shows the difference between the input and output angles 
versus input angle for the HS-A incrossing angle. 
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Fig. 2. The OrhView-2 horizon scanner geometry. 


Each angle is then converted to a unit horizon crossing 
vector in the spacecraft frame: 


and 


Sh 


cos $ sin 7 
sin ^ sin 7 , 
cos 7 


Us = |H| S H , 


(42) 


(43) 


where 7 is the scanner half-cone angle of 45° , Hh is the 
horizon vector in the scanner reference frame, Hs is the 
vector in the spacecraft frame, and H is the scanner frame- 
to-spacecraft frame transformation matrix. 


3 .2. 2. 2 Horizon Angle Model 

The model horizon height at each crossing is a function 
of several variables; the orbit geocentric latitude and posi- 
tion magnitude, the azimuth of the crossing point, and the 
height of the C0 2 layer (Liu 1978). Standard equations are 
available for computing the horizon angle based on static 
models; however, analysis of the OV2 horizon scanner data 
has shown that the apparent CO 2 height shifts seasonally 
(Patt and Bilanow 2001). The model described below was 
derived for SeaWiFS navigation in order to allow for this 
shift. 

The horizon angle to the geocentric nadir, p, is com- 
puted by solving the following quadratic equation: 

0 = iJ 2 cos 2 ^ - 4iT 6 sin 2l I r - 4H x He cos 2 $ 

+ (4ff 3 tf 6 cos - 2H4H5 cos $) cot p (44) 
+ {Hi - 4J7 2 tf 6 )cotV 


where $ is the azimuth of the horizon crossing relative to 
local north, and the coefficients H x through H 6 are com- 


puted as follows: 



Hi = 

2 cos 2 A 

Sm A+ (l_/)2’ 


(45) 

H a = 

2 sin 2 A 

COS A + ^ ’ 


(46) 


. . 2 sin A cos A 


(47) 

H 3 = 

—2 sin A cos A H 7- , 

(1 -Jr 



. - 21^1 sin A cos A - 2Szcos A 

,(48) 


— 2|P| sin A cos A -I — 

~f ) 2 

H 5 = 

- « 2\P\ sin 2 A — 2Sz sin A 

2|P| cos 2 A + ‘ ' (1 _ /)2 

(49) 

and 




#6 = 


(Re + Sr) 2 , 

(50) 


where A is the orbit geocentric latitude, Sz is a seasonal 
ellipsoid shift along the Earth’s pole, and Sr is a seasonal 
adjustment to the Earth’s radius. The value of / for the 
Earth’s surface is 1/298.257 (Wertz 1978); however, the 
flattening of the C0 2 layer is adjusted based on the anal- 
ysis of the horizon scanner data (Patt and Bilanow 2001): 


1 sin A s 
185 64Cf’ 


(51) 


where As is the subsolar latitude for the sample time. The 
remaining seasonal corrections are computed as: 


Sz = 11 sin As, (52) 

and 
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Fig. 3. The in-crossing angle calibration for HS-A, where the units are in degrees. 


Sr = — 4. 5 sin As. 


(53) 


The horizon crossing azimuth, T, depends on the az- 
imuth of the orbit reference frame y- axis, the spacecraft 
yaw angle, the scanner rotation axis in the spacecraft frame, 
and the Earth chord measurement. The orbital y-axis az- 
imuth is determined relative to the local north direction at 
the orbit position, which is computed as 


ft = 


' -P X Pz ■ 
-PyP* 
Pl + Py 


1 

yJ\P\ (Pf+Pf) 


(54) 


where $ is the unit vector in the northern direction. The 
orbital azimuth given by 

(55) 

where To is the azimuth of the orbital frame y-axis, and 
O y and O z are the orbital y- and z-axes, respectively. The 
spacecraft y-axis azimuth differs from the orbited azimuth 
by the yaw angle, and the scanner axis azimuths are ±90° 
from the spacecraft, so the scanner axis azimuth is given 

by 

Tj/s = T o +i/>-90, (56) 


for HS-A and the same equation with +90° is used for HS- 
B. The next step is to compute the azimuth offset from the 
scanner rotation axis to the horizon crossings. This offset 
is computed from the chord by first computing the horizon 
angle p 0 at the scanner axis azimuth using (45). Using this 
angle, the chord, and the scanner cone angle are used to 
compute the azimuth offset: 

. _ _i I sin : 

AT = cos < — 


7 cos ft + cos 2 7 — cos 2 p 
2 sin 2 p 


(57) 


and 

T = + AT, (58) 


where the offset is positive for the in-crossing and the same 
equation is used with a negative offset for the out-crossing. 

The final step is to adjust p to account for the CO2 
layer height. This height is assumed to be 40 km, which 
is included by adding a bias of 0.74° to the value of p 
determined from (44). 

Note that the spacecraft yaw angle is needed for (56). 
This would seem to be circular logic, because the nadir 
vector is ultimately used to compute yaw. The horizon 
angle, however, varies slowly with yaw (about 0.003), so 
the input to (56) need only be accurate to a few degrees. 
For the first scan, the yaw is assumed to be zero, and the 
calculations are repeated if the final yaw is larger than a 
specified tolerance. For all subsequent scans, the yaw is 
used from the previous scan. 
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3. 2. 2. 3 Nadir Vector Calculation 

The calculation of the nadir vector can be performed 
using data from either one or both scanners. If one scanner 
is used, the calculation is deterministic; for two scanners 
a least-squares estimation is performed. In either case the 
calculation is based on the following relationship between 
a horizon crossing vector, the corresponding horizon angle 
and the nadir vector: 


H E = cos /9, 


and 


i?l + &2 

I#1 + s 2 f 

= fflX S 2 
2 l&x S 2 f 

Si - h 2 


Ui = 


tJ 2 = 


& = 


\Si - S 2 \ 


(60) 


(61) 


(62) 


Combining (60) and (62) with (59) above produces the 
following: 


jj g = COS Pi + COSP2 

1 ’ \ffi + S 2 \ 


= u 1, 


and 


tJ 3 ■£ = 


COS P\ 


COS P 2 


\H X - H 2 1 


U 3- 


(63) 


(64) 


By definition, because t) i, and tJ 3 are orthonormal, 


)/l - u l - u l 

U 2 , 

(65) 

U\0\ + U 2 O 2 "h U3^?3, 

(66) 


and, therefore, 


where 14, U2, and 1x3, are intermediate coefficients used to 
calculate the nadir vector, and the sign of U2, is chosen to 
ensure that the x component of the nadir vector is positive. 

For more than two horizon crossing vectors (dual- 
scanner processing), the nadir vector is over determined, 
and a least-squares calculation is performed. The loss func- 
tion is constructed from (59) as follows: 


J = 0.5V [ff.E — COS Pi] 2 . 


(67) 


Minimizing (67) with respect to E yields: 

4 4 

£ 

1 = 1 


(59) 


where £ is the geocentric nadir vector in the spacecraft 
frame. 

The deterministic calculation is performed by con- 
structing a set of orthonormal axes using the two horizon 
crossing vectors, and then applying the above relationship 
to solve for E. The axes are constructed as follows: 


cos Pi =0. (68) 

1=1 t=l 

The formulation in (68) can be rearranged: 

cos p*, (69) 

which is solved for & by matrix inversion and multiplica- 
tion. The final nadir vector is normalized to correct for 
any numerical errors. 

4. ATTITUDE DETERMINATION 

The attitude determination processing for SeaWiFS 
uses the computed sun and nadir vectors in the body frame, 
along with the reference sun and nadir vectors in the ECEF 
frame, to compute the yaw, roll, and pitch angles. A 
Kalman smoothing algorithm is used to reduce the un- 
certainty in the yaw angle near the subsolar point, allow 
processing through gaps in attitude sensor data, and min- 
imize discontinuities caused by transitions in sensor cov- 
erage (i.e., DSS FOVs). The dynamics model is based on 
the assumption that the pitch axis is nearly inertial, which 
is due to the pitch momentum bias resulting from the mo- 
mentum wheel. 

The following sections present the equations used, fol- 
lowed by a description of the implementation. 

4.1 Kalman Filter Equations 

The basic Kalman filter equations are given (Fallon 
1978) as follows. The state update equations are: 


X(t) = X-i (t) + 


G), 


= P_i(*)G t [R + GP_i(£)G 


T\~ 1 


(70) 

(71) 


and 


P = [n- 


_i[l-KG] r + KRK r , (72) 


where X (t) is the state vector of dimension 3 at time t\ Y is 
the observation vector of dimension 6; G is the observation 
model vector, also of dimension 6; K is the (3 x 6) gain 
matrix; P is the (3 x 3) state covariance matrix; G is the 
(6 x 3) matrix of partial derivatives of G with respect to 
X] and R is the (6 x 6) observation covariance matrix. 
The subscript (—1) on X and P indicates a state at time t 
propagated from the previous time. 

The state propagation equations are: 




X(t + At) = D (t + At,t)X(t), 


(73) 
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and 

P(t + At) = B(t + At,*) P(t) 0 T (f + At, t) 

+ Q(t + At, t), 

where D is the (3x3) state transition matrix, Q is the 
(3 x 3) state noise covariance matrix, and t is the time 
step. 

For SeaWiFS attitude determination, the state X con- 
sists of the yaw, roll, and pitch attitude angles, and the 
observations V consist of the components of the computed 
sun and nadir vectors. The calculations of G , G, R, Q, 
and P are performed as described in the remainder of this 
section. 

The observation model is based on the reference sun 
and nadir vectors transformed to the spacecraft frame. 
The sun reference vector, Se , is computed in the ECEF 
frame from a model (Patt and Gregg 1994) and trans- 
formed from ECEF-to-spacecraft coordinates using the 
transformations given in Sect. 1.2: 


-sin# — cos# 0 

— = cos # —sin# 0 . (81) 

09 [ 0 0 0_ 

The 6x6 observation covariance matrix R is specified as 
diagonal, a simplification which ignores any coupling be- 
tween the vector components. The covariances are based 
on an estimated uncertainty of 0.1° for the horizon scan- 
ners and 0.06° for the sun sensors. The sun sensor un- 
certainties are assumed to be constant in angular terms, 
which means the vector component variances decrease as 
the components (S x , S yy and S z ) approach dbl. The matrix 
is specified as: 

R Sx 0 0 0 0 0 

0 R Sy 0 0 0 0 

0 0 R s , 0 0 0 

0 0 0 Rh 0 0 

0 0 0 0 R h 0 

0 0 0 0 0 Rh 




Sr = MOSe, (75) 

where Sr is the sun reference vector in the spacecraft 
frame; and B is computed using the propagated yaw, roll, 
and pitch. The reference nadir vector is the negative of the 
normalized orbit position vector, which is also transformed 
to spacecraft coordinates: 



and, therefore, 



(77) 


The calculations for determining G, the partial derivatives 
of these vectors with respect to can be performed by 
recognizing that only the orbital-to-spacecraft transforma- 
tion depends on X\ e.g., 


OSr 

dX 



(78) 


where Rs x , Rs y , and Rs , are the variances in S x , S y , and 
S Z) respectively, and Rh is the variance in the horizon 
vector, based on the above-stated uncertainties. 

The 3x3 state noise covariance matrix Q is also diag- 
onal, and is based on the expected variance of the actual 
attitude from the simple dynamics model. It is expressed 
as; 

q^At 0 0 

Q = 0 q<j>At 0 , (83) 

0 0 q&At_ 

where q<p, and q& are the estimated time variances in 
yaw, roll, and pitch, respectively, and At is the time be- 
tween scan lines, i.e. , 2/3 s for GAC data and l/6s for 
LAC and HRPT. The variances need to allow for all pitch 
motion, which is not modeled, and also for the deviation 
of roll and yaw from the model of the inertial spin axis. 
The values were all chosen, after some experimentation, to 
be 2.5 x 10 -9 s _1 . 

The state transition matrix, D, is based on the dynam- 
ics model of the inertial spin axis, i.e., constant pitch angle 
and orbital coupling between roll and yaw. This is speci- 
fied using a small angle approximation as 


and likewise for Sr. It is only necessary, therefore, to 
differentiate the expression for B, (8), with respect to the 
yaw, roll, and pitch angles, as follows: 


D = 


1 -w 0 At 0" 
bj 0 At 1 0 , 


(84) 


and 



0 0 

0 


0 — sin ip 

cosip A^Afl, 

(79) 

0 — cos rp 

— sin ip 


— sin 0 0 

cos <p 


0 0 

0 As, 

(80) 

— cos (p 0 

— sin <f> 



4.2 Kalman Smoother Implementation 

The Kalman filter, as described by Fallon (1978), uses 
only previous data to estimate the state. The SeaWiFS 
attitude determination processing can take advantage of 
the fact that, except at the ends of a scene, both past and 
future measurements are available. This is done by making 
two passes through the scene, forwards and backwards, and 
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then computing a weighted average of the results for each 
scan line — this process is known as Kalman smoothing. 

The processing is initialized for a scene by setting the 
state vector, A*, to all zeroes and P, the state covariance 
matrix, to an identity matrix. For each scan line, the pro- 
cessing is as follows. 

1. The state transition matrix and state noise covari- 
ance matrix (B and Q, respectively) are computed 
using (83) and (84) based on the time from the pre- 
vious scan line, At. 

2. The previous values of X and P are then propagated 
using (73) and (74). 

3. Then, if valid observations of Y are available, the 
model observations, partial derivatives, and obser- 
vation covariance matrix (<5, G, and M, respectively) 
are computed using (75) through (82). 

4. The gain matrix K is computed using (71), and is 
used to update X and P using (70) and (72). 

5. The values of X and P are saved for each line. 

The second pass is performed using X and P from the 
last scan line for initialization. The computations are iden- 
tical except that At is negative. At each line, the weighted 
average of X is computed using the values of P from both 
passes, and saved as the final set of spacecraft attitude 
angles for that scan line. 

5. FINAL PROCESSING 

The final steps in navigation processing are: processing 
the SeaWiFS instrument tilt telemetry to determine the 
tilt angle for each scan line; combining the orbit, attitude 
and tilt information to compute the instrument pointing 
information for each scan line; and generating the data to 
be stored in the level -la output product. These steps are 
described in the following subsections. 

5.1 Tilt Data Processing 

The SeaWiFS instrument can be commanded to three 
tilt states: nadir (0°), aft (20°), and forward (—20°). As 
stated in Sect. 1, the tilt consists of a rotation about the 
SeaWiFS z-axis. The standard operating scenario is to 
tilt aft for the first half of the data collection period in the 
orbit (i.e., north of the subsolar point) and forward for the 
second half. Other scenarios, primarily involving nadir tilt 
for part of an orbit, have been used on occasion. The tilt 
change takes about 13 s to complete. 

The tilt telemetry consists of digital status bits and 
analog tilt motor angles, which are provided for two out of 
every three scan lines. The three status bits indicate one of 
the four possible states: nadir, aft, forward, or changing. 
In the first three states, the tilt angle can be assumed to be 
a fixed, predetermined value. The tilt angles at the three 
positions were measured prelaunch and is given in Table 2 


(note that nadir has been defined to be exactly zero). Dur- 
ing the tilt change periods, the tilt motor telemetry is used 
to compute the tilt angle, as described below. The tilt 
telemetry are validated by checking for single scan-line tilt 
states (indicative of corrupted telemetry), out-of-limit val- 
ues, and monotonic behavior during tilt states. The tilt 
angles are interpolated to the scan lines without tilt infor- 
mation for which the tilt telemetry were flagged. 


Table 2. Nominal and actual tilt angles 


Tilt State 

Nominal Angle 

Actual Angle 

Nadir 

0° 

o 

O 

O 

Aft 

o 

O 

CN 

19.820° 

Forward 

-20° 

-19.925° 


The tilt motor telemetry consists of angle measure- 
ments from two tilt motors, one mounted on the fixed base 
of the sensor and one on the tilting platform. The angles 
are referred to as and $p, respectively. After valida- 
tion, the next step is to calibrate the angles and add an 
adjustment to the measurement reference points: 


= Ab$b + B b 4 - $RBy ( 85 ) 


and 

= Ap$p + Bp + 4>pp, (86) 

where and $p are the calibrated and adjusted motor 
angles; A b and Ap are the calibration scale factors; B b , 
and Bp are the biases; and <&r B and $pp are the reference 
angle adjustments (80 and 235°, respectively). 

The tilt mechanism includes the following fixed dimen- 
sions: 


a) Pivot- to-base motor shaft, d B = 7.088; 

b) Pivot-to-platform motor shaft, dp = 6.000; 

c) Base-to-platform motor link, di = 4.070; and 

d) Motor shaft- to-link radius, d r = 1.025. 

The tilt angle, 0, is computed as follows: 

0 =©' - ©p - ©p - 35.068, (87) 


where 


©' 

©B 

©p 


cos 


' d' B 2 + d'p 2 - df 


sin 


(2 d'gd'p) 
_1 / dr sin \ 


sin 


-l 


d r sin $>p 

dX, 


(88) 

(89) 

(90) 


d'g = \Jd? B + ~d% - 2d B d T cos$' B , (91) 

and 

dp = yjd 2 p + d 2 — 2dp d r cos 4>p, (92) 
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5.2 Final Attitude Processing 

The data stored for each scan line consist of the orbit 
position vector, attitude angles, tilt angle, sun reference 
vector in the ECEF frame, sensor (ECEF- to- SeaWiFS) 
transformation matrix, and scan ellipse coefficients. In ad- 
dition, flags are set if navigation cannot be computed for 
a scan line. 

The computation of the orbit vector, attitude angles, 
tilt angle, and sun reference vector were described in pre- 
vious sections. The sensor transformation matrix is com- 
puted by successively multiplying the ECEF-to-orbital, 
orbital-to-spacecraft, spacecraft-to-SeaWiFS, and tilt 
transformations: 

M = TSIO, (93) 

where the orbital-to-spacecraft transformation is computed 
using (8) and the yaw, roll, and pitch angles for the scan 
line. 

The six scan ellipse coefficients are computed using 
the sensor transformation matrix and orbit position vector 
(Patt and Gregg 1994): 


Ci 


C 2 


C 3 

C 4 

C 5 


+ M ?,» + (fT/)5 
i + m;, 3 ((i - /)- 2 - i), 


= + 2Mi ; 2M3 i 2 + 


2Ml ; 3M; 


3,3 


(i - /) 2 


= 2M 1i3 M 3 ,3 ((1 - /)- 


1 ). 


M§,3 


M 3 ,i + M 3> 2 + ^ _ ' jy 

1 + m 3 3 ((1 - /)- 2 - 1), 

2M 1i3 P, 


= 2Mi i iPr + + 

= 2M 3> iPc + 2M 3j 2Py + 


(i - tr 
2M 3i3 P, — 
(1 - /) 2 ’ 


(94) 

(95) 

(96) 

(97) 

(98) 


and 

C 6 = Pi + Pi + P ' /)2 - R\, (99) 

where Mi t j is the component of M in the ith row and jth 
column. 


6. DISCUSSION 

The navigation algorithms presented here are imple- 
mented in the operational data processing code and are 
used daily in the data processing performed by the Sea- 
WiFS Projectf. The accuracy of the navigation accu- 
racy has been assessed continually through an automated 


method using island targets (Patt et al. 1997). This assess- 
ment has served to provide a public record of the accuracy 
of the current navigation processing!, to indicate where 
the accuracy needs improvement, and to assess the effect 
of changes to the algorithms prior to their acceptance for 
operational processing. 

Provided below is a brief summary of the accuracy as- 
sessed as various improvements to the processing were in- 
troduced since launch. It includes the typical accuracies 
now expected of the algorithms, and some cases that pro- 
vide exceptions to the nominal performance are noted. In 
addition, brief observations are made on the potential for 
further improvements. 

The changes to the algorithms since launch produced 
rapid improvements in accuracy early in the mission, and 
smaller improvements over longer time scales since then. 
Shortly after launch, the initial adjustments to the process- 
ing resulted in a typical 2a accuracy of about 2 pixels, with 
considerable variation; the errors near the subsolar point 
were largely unknown. In early 1998, the attitude sensor 
alignments were refined, which reduced the typical accu- 
racy to about 1.5 pixels, but the subsolar errors remained 
large. 

During mid-1998, the Kalman smoother was imple- 
mented for attitude determination to replace the single- 
frame method, which was developed before launch. This 
was performed in anticipation of the single- ACS-string op- 
eration, which was planned to start later that year, and 
also to improve the subsolar accuracy. In this same time 
frame, the horizon angle model was modified to include 
a revised flattening factor and seasonal corrections. The 
combination of these improvements and additional refine- 
ment of the attitude sensor alignments produced typical 
accuracies of about 1.25 pixels, with subsolar errors of up 
to several pixels. These changes were incorporated for the 
reprocessing of the mission data in August and September 
1998. 

In 1999, analysis of the navigation assessment results 
for the entire mission indicated inconsistencies in the roll 
and yaw accuracies, with better results during the months 
of April-September at the expense of larger errors during 
October-March. This was traced to the use of a simpli- 
fied form of the horizon scanner angle calibration, which 
was implemented shortly after launch. The calibration was 
modified to use look-up tables provided by the manufac- 
turer. This improvement was used for the reprocessing in 
May 2000. With this change, the typical accuracies ap- 
proached 1 pixel; errors of a few pixels were seen near, 
and just south of, the subsolar point, particularly during 
December-March . 

Most recently, in early 2001, a further refinement to 
the horizon angle model was developed in response to nav- 
igation assessment results which indicated seasonal varia- 
tions in the horizon flattening factor. A correction for this 


f The SeaWiFS processing code is freely distributed by NASA 
to SeaWiFS data users and ground stations as part of the 
SeaWiFS Data Analysis System (SeaDAS). 


f Available via the World Wide Web at the following URL: 
http : \\algae2 . gsf c . nasa . gov/navqc . 
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was implemented, and the navigation accuracies since then 
have averaged 1 pixel or less, with maximum errors of less 
than 2 pixels. This change will be applied to the mission 
data set during the reprocessing in March or April 2002. 

The typical accuracy is not achieved in a variety of 
special cases: 

• Time tag errors; Occasional errors in the time tags 
for the science data have been discovered, which 
offset the geolocation in the along-track direction. 
Efforts are underway to correct or flag these occur- 
rences during level- la processing. 

• Polar horizon anomalies: Winter polar stratospheric 
temperature variability affects the horizon scanner 
performance at high latitudes, degrading the roll 
accuracy. 

• Periods of high nutation: High levels of nutation 
were triggered which occurred often in the early 
part of the mission, and occasionally thereafter. At 
these times, the tracking of the roll and yaw changes 
was less accurate. 

• Anomalous pitch motion: Some fairly rapid pitch 
motions of the spacecraft have been caused by on- 
board GPS errors and from moon interference in the 
scanners; the current processing algorithms cannot 
always track these motions well. 

• Tilt change: During the SeaWiFS instrument tilt 
change and for a few seconds before and after, the 
pitch accuracy is affected. 

• Short data spans: For HRPT scenes less than about 
5 min duration, especially in the subsolar region, ac- 
curacy is degraded. 

Although further improvements are certainly possible, 
they would mainly reduce the maximum errors under cer- 
tain unusual conditions as discussed above. The accuracy 
is currently pushed near the limits set by the performance 
of the attitude sensors and the fidelity of the spacecraft 
dynamics model. Further improvements would probably 
need to depend on modeling the spacecraft dynamics and 
onboard control responses, and this would present some 
challenges. Ultimately, this activity will continue accord- 
ing to the needs of the science-data users. 


Glossary 

ACS Attitude Control System 
ASAP Artificial Satellite Analysis Program 

COSMIC Computer Software Management and Information 
Center 

DSS Digital Sun Sensor 

ECEF Earth- Centered Earth-Fixed 
ECI Earth- Centered Inertial 
ECR Earth- Centered Rotating 

FOV Field-Of-View 


GAC Global Area Coverage 
GPS Global Positioning System 
GSFC Goddard Space Flight Center 

HRPT High-Resolution Picture Transmission 
HS Horizon Scanner 

IQR Interquartile Range 

LAC Local Area Coverage 

NASA National Aeronautics and Space Administration 

OV 2 OrbView- 2 ; satellite platform on which the Sea- 
WiFS instrument is flown. 

SeaDAS SeaWiFS Data Analysis System 
SeaWiFS Sea- viewing Wide Field-of-view Sensor 

URL Universal Resource Locator 
Symbols 

a Orbit semi-major axis. 

Ab Tilt motor angle scale factor. 

Ap Tilt motor angle scale factor. 

Acc A DSS scale factor. 

Ap A DSS scale factor. 

A Roll angle rotation matrix. 

Pitch angle rotation matrix. 

Av> Yaw angle rotation matrix. 

B Orbital- to-spacecraft transformation matrix. 

Bq Tilt motor angle bias. 

Bp Tilt motor angle bias. 

B a DSS bias. 

Bp DSS bias, 

cb Cubic polynomial coefficient vector, 
ci Cubic polynomial coefficient vector. 

C2 Cubic polynomial coefficient vector. 

C3 Cubic polynomial coefficient vector. 

C\ Scan ellipse coefficient. 

C2 Scan ellipse coefficient. 

C3 Scan ellipse coefficient. 

C4 Scan ellipse coefficient. 

C5 Scan ellipse coefficient. 

Ce Scan ellipse coefficient, 

C DSS- to-spacecraft transformation matrix. 

D Kalman filter state transition matrix. 
ds Tilt pivot to base motor shaft dimension. 
d' B Tilt pivot to base motor link dimension. 
d\ Tilt base motor to platform motor link dimension. 
dp Tilt pivot to platform motor shaft dimension. 
d f p Tilt pivot to platform motor link dimension. 
d r Tilt motor shaft to link dimension, 

e Orbit eccentricity. 

S Nadir vector in spacecraft frame. 

E ECI-to-ECEF transformation matrix. 

Sr Reference nadir vector. 

/ Earth flattening factor. 

fp Effective flattening factor at spacecraft position. 

G Observation model vector. 

G Partial derivative matrix. 

G m Earth gravitational constant. 
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ff Horizon vector in spacecraft frame. 

EH HS- to- spacecraft transformation matrix. 

Hi Horizon angle model coefficient. 

H 2 Horizon angle model coefficient. 

H3 Horizon angle model coefficient. 

H4 Horizon angle model coefficient. 

H5 Horizon angle model coefficient. 

He Horizon angle model coefficient. 

I$h Horizon vector in HS frame. 

Us Horizon vector in spacecraft frame. 

i Orbit inclination. 

I Identity matrix. 

J Loss function. 

J 2 Earth gravity field perturbation term. 

K Kalman filter gain matrix. 

I Orbit right ascension of ascending node. 

L Orbit plane normal vector. 

L x x component with relation to L. 

L y y component with relation to L, 

L z z component with relation to L. 
m Orbit mean anomaly. 

M ECEF-to-SeaWiFS frame transformation matrix, 
n Mean motion. 

$ North vector at orbit position. 
o Orbit angle measured from the ascending node. 
O ECEF-to-orbital transformation matrix. 

O x Orbital frame axis. 

O v Orbital frame axis. 

O z Orbital frame axis. 

P Orbit position vector. 

Pi Position vector at time t\. 

P 2 Position vector at time t 2 . 

P State covariance matrix. 

P x x component with relation to P. 

P y y component with relation to P. 

P z z component with relation to P. 

qe Pitch angle state noise covariance. 
q<f, Roll angle state noise covariance. 
q Yaw angle state noise covariance. 

Q State noise covariance matrix. 

M Kalman filter measurement covariance matrix. 
Re Earth equatorial radius. 

Rh Horizon scanner measurement variance. 

Rm Earth mean radius. 

Rs Sun sensor measurement variance. 

§ Spacecraft-to-SeaWiFS transformation matrix. 
Pd Sun vector in DSS frame. 

Se Sun reference vector in ECEF frame. 

P R Sun reference vector in spacecraft frame. 

Ps Sun vector in spacecraft frame^ 

S x x component with relation to S. 

S y y component with relation to S. 

S z z component with relation to S. 
t Relative sample time. 
t s Sample time, 
to Epoch time. 

t\ Time of first vector used for interpolation. 
t 2 Time of second vector used for interpolation. 

T Transpose (of a matrix). 

T Tilt transformation matrix. 


u\ Intermediate coefficient used to calculate the nadir 
vector. 

U 2 Intermediate coefficient used to calculate the nadir 
vector. 

U 3 Intermediate coefficient used to calculate the nadir 
vector. 

Vi Intermediate vector for nadir vector calculation. 

02 Intermediate vector for nadir vector calculation. 

03 Intermediate vector for nadir vector calculation. 

V Orbit velocity vector. 

Pi Velocity vector at time t \ . 

V 2 Velocity vector at time £ 2 . 

Vc Orbit velocity vector in ECEF frame corrected for 
Earth rotation rate. 

V x x component with relation to V'. 

V y y component with relation to 

V x z component with relation to P. 

w Orbit argument of perigee. 

W State weighting matrix. 

x Reference frame axis, or component of a vector in 
this frame. 

X State vector. 

X(t) State vector at time t. 

y Reference frame axis, or component of a vector in 
this frame. 

P Observation vector (GPS orbit position or ACS sen- 
sor vector measurements) . 

z Reference frame axis, or component of a vector in 
this frame. 

P Geocentric North Pole vector. 

a DSS output angle. 

cl DSS calibrated output angle. 

P DSS output angle. 

P f DSS calibrated output angle. 

7 HS half-cone angle. 

Sr Seasonal adjustment to the Earth’s radius. 

Sz Seasonal ellipsoid shift along the Earth’s pole. 

At Time difference. 

AX State update vector. 

Azimuth offset. 

6 Spacecraft pitch angle. 

0 SeaWiFS tilt angle. 

0' Intermediate component of tilt calculation from ( 88 ) 
to (90). 

0b Intermediate component of tilt calculation from ( 88 ) 
to (90). 

Op Intermediate component of tilt calculation from ( 88 ) 
to (90). 

A Subsatellite latitude. 

A s Subsolar latitude. 

a Standard deviation. 

v Orbit true anomaly. 

£ Greenwich Hour Angle. 

<f> Spacecraft roll angle. 

3>b Raw tilt motor angle. 

$>' B Calibrated tilt motor angle. 

HS crossing angle (in-crossing). 

$0 HS crossing angle (out-crossing). 

Raw tilt motor angle. 
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$p Calibrated tilt motor angle. 

Tilt motor reference angle. 

Tilt motor reference angle. 

4>s HS phase angle output. 
p Horizon angle to geocentric nadir. 
ip Spacecraft yaw angle. 

^ Azimuth of the horizon crossing relative to local 
North. 

$hs Azimuth of HS axis. 

Azimuth of orbital y axis. 

uje Earth rotation rate. 
ujo Orbit rotation rate. 

Q Horizon scanner Earth chord angle output. 
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